High-Throughput Sequencing of microRNAs in Peripheral Blood Mononuclear Cells: Identification of Potential Weight Loss Biomarkers

INTRODUCTION: MicroRNAs (miRNAs) are being increasingly studied in relation to energy metabolism and body composition homeostasis. Indeed, the quantitative analysis of miRNAs expression in different adiposity conditions may contribute to understand the intimate mechanisms participating in body weight control and to find new biomarkers with diagnostic or prognostic value in obesity management. OBJECTIVE: The aim of this study was the search for miRNAs in blood cells whose expression could be used as prognostic biomarkers of weight loss. METHODS: Ten Caucasian obese women were selected among the participants in a weight-loss trial that consisted in following an energy-restricted treatment. Weight loss was considered unsuccessful when <5% of initial body weight (non-responders) and successful when >5% (responders). At baseline, total miRNA isolated from peripheral blood mononuclear cells (PBMC) was sequenced with SOLiD v4. The miRNA sequencing data were validated by RT-PCR. RESULTS: Differential baseline expression of several miRNAs was found between responders and non-responders. Two miRNAs were up-regulated in the non-responder group (mir-935 and mir-4772) and three others were down-regulated (mir-223, mir-224 and mir-376b). Both mir-935 and mir-4772 showed relevant associations with the magnitude of weight loss, although the expression of other transcripts (mir-874, mir-199b, mir-766, mir-589 and mir-148b) also correlated with weight loss. CONCLUSIONS: This research addresses the use of high-throughput sequencing technologies in the search for miRNA expression biomarkers in obesity, by determining the miRNA transcriptome of PBMC. Basal expression of different miRNAs, particularly mir-935 and mir-4772, could be prognostic biomarkers and may forecast the response to a hypocaloric diet.

SNP Discovery in European Anchovy (Engraulis encrasicolus, L) by High-Throughput Transcriptome and Genome Sequencing

12 p.

A short in-frame deletion in NTRK1 tyrosine kinase domain caused by a novel splice site mutation in a patient with congenital insensitivity to pain with anhidrosis

Background: Congenital insensitivity to pain with anhidrosis (CIPA) is a rare autosomal recessive genetic disease characterized by the lack of reaction to noxious stimuli and anhidrosis. It is caused by mutations in the NTRK1 gene, which encodes the high affinity tyrosine kinase receptor I for Neurotrophic Growth Factor (NGF). -- Case Presentation: We present the case of a female patient diagnosed with CIPA at the age of 8 months. The patient is currently 6 years old and her psychomotor development conforms to her age (RMN, SPECT and psychological study are in the range of normality). PCR amplification of DNA, followed by direct sequencing, was used to investigate the presence of NTRK1 gene mutations. Reverse transcriptase (RT)-PCR amplification of RNA, followed by cloning and sequencing of isolated RT-PCR products was used to characterize the effect of the mutations on NTRK1 mRNA splicing. The clinical diagnosis of CIPA was confirmed by the detection of two splice-site mutations in NTRK1, revealing that the patient was a compound heterozygote at this gene. One of these alterations, c.574+1G > A, is located at the splice donor site of intron 5. We also found a second mutation, c.2206-2 A > G, not previously reported in the literature, which is located at the splice acceptor site of intron 16. Each parent was confirmed to be a carrier for one of the mutations by DNA sequencing analysis. It has been proposed that the c.574+1G > A mutation would cause exon 5 skipping during NTRK1 mRNA splicing. We could confirm this prediction and, more importantly, we provide evidence that the novel c.2206-2A > G mutation also disrupts normal NTRK1 splicing, leading to the use of an alternative splice acceptor site within exon 17. As a consequence, this mutation would result in the production of a mutant NTRK1 protein with a seven aminoacid in-frame deletion in its tyrosine kinase domain. --Conclusions: We present the first description of a CIPA-associated NTRK1 mutation causing a short interstitial deletion in the tyrosine kinase domain of the receptor. The possible phenotypical implications of this mutation are discussed.

Genetics of Type III Bartter Syndrome in Spain, Proposed Diagnostic Algorithm

9 p.

Molecular diagnosis of distal renal tubular acidosis in Tunisian patients: proposed algorithm for Northern Africa populations for the ATP6V1B1, ATP6V0A4 and SCL4A1 genes

Background: Primary distal renal tubular acidosis (dRTA) caused by mutations in the genes that codify for the H+ -ATPase pump subunits is a heterogeneous disease with a poor phenotype-genotype correlation. Up to now, large cohorts of dRTA Tunisian patients have not been analyzed, and molecular defects may differ from those described in other ethnicities. We aim to identify molecular defects present in the ATP6V1B1, ATP6V0A4 and SLC4A1 genes in a Tunisian cohort, according to the following algorithm: first, ATP6V1B1 gene analysis in dRTA patients with sensorineural hearing loss (SNHL) or unknown hearing status. Afterwards, ATP6V0A4 gene study in dRTA patients with normal hearing, and in those without any structural mutation in the ATP6V1B1 gene despite presenting SNHL. Finally, analysis of the SLC4A1 gene in those patients with a negative result for the previous studies. Methods: 25 children (19 boys) with dRTA from 20 families of Tunisian origin were studied. DNAs were extracted by the standard phenol/chloroform method. Molecular analysis was performed by PCR amplification and direct sequencing. Results: In the index cases, ATP6V1B1 gene screening resulted in a mutation detection rate of 81.25%, which increased up to 95% after ATP6V0A4 gene analysis. Three ATP6V1B1 mutations were observed: one frameshift mutation (c.1155dupC; p.Ile386fs), in exon 12; a G to C single nucleotide substitution, on the acceptor splicing site (c.175-1G > C; p.?) in intron 2, and one novel missense mutation (c. 1102G > A; p. Glu368Lys), in exon 11. We also report four mutations in the ATP6V0A4 gene: one single nucleotide deletion in exon 13 (c.1221delG; p. Met408Cysfs* 10); the nonsense c.16C > T; p.Arg6*, in exon 3; and the missense changes c.1739 T > C; p.Met580Thr, in exon 17 and c.2035G > T; p.Asp679Tyr, in exon 19. Conclusion: Molecular diagnosis of ATP6V1B1 and ATP6V0A4 genes was performed in a large Tunisian cohort with dRTA. We identified three different ATP6V1B1 and four different ATP6V0A4 mutations in 25 Tunisian children. One of them, c.1102G > A; p.Glu368Lys in the ATP6V1B1 gene, had not previously been described. Among deaf since childhood patients, 75% had the ATP6V1B1 gene c. 1155dupC mutation in homozygosis. Based on the results, we propose a new diagnostic strategy to facilitate the genetic testing in North Africans with dRTA and SNHL.

Multi-drug resistance profiles and the genetic features of Acinetobacter baumannii isolates from Bolivia

Introduction: Acinetobacter baumannii is opportunistic in debilitated hospitalised patients. Because information from some South American countries was previously lacking, this study examined the emergence of multi-resistant A. baumannii in three hospitals in Cochabamba, Bolivia, from 2008 to 2009. Methodology: Multiplex PCR was used to identify the main resistance genes in 15 multi-resistant A. baumannii isolates. RT-PCR was used to measure gene expression. The genetic environment of these genes was also analysed by PCR amplification and sequencing. Minimum inhibitory concentrations were determined for key antibiotics and some were determined in the presence of an efflux pump inhibitor, 1-(1-napthylmethyl) piperazine. Results: Fourteen strains were found to be multi-resistant. Each strain was found to have the bla(OXA-58) gene with the ISAba3-like element upstream, responsible for over-expression of the latter and subsequent carbapenem resistance. Similarly, ISAba1, upstream of the bla(ADC) gene caused over-expression of the latter and cephalosporin resistance; mutations in the gyrA(Ser83 to Leu) and parC (Ser-80 to Phe) genes were commensurate with fluoroquinolone resistance. In addition, the adeA, adeB efflux genes were over-expressed. All 15 isolates were positive for at least two aminoglycoside resistance genes. Conclusion: This is one of the first reports analyzing the multi-drug resistance profile of A. baumannii strains isolated in Bolivia and shows that the over-expression of thebla(OXA-58), bla(ADC) and efflux genes together with aminoglycoside modifying enzymes and mutations in DNA topoisomerases are responsible for the multi-resistance of the bacteria and the subsequent difficulty in treating infections caused by them.

Assembly of viral genomes from metagenomes

Viral infections remain a serious global health issue. Metagenomic approaches are increasingly used in the detection of novel viral pathogens but also to generate complete genomes of uncultivated viruses. In silico identification of complete viral genomes from sequence data would allow rapid phylogenetic characterization of these new viruses. Often, however, complete viral genomes are not recovered, but rather several distinct contigs derived from a single entity are, some of which have no sequence homology to any known proteins. De novo assembly of single viruses from a metagenome is challenging, not only because of the lack of a reference genome, but also because of intrapopulation variation and uneven or insufficient coverage. Here we explored different assembly algorithms, remote homology searches, genome-specific sequence motifs, k-mer frequency ranking, and coverage profile binning to detect and obtain viral target genomes from metagenomes. All methods were tested on 454-generated sequencing datasets containing three recently described RNA viruses with a relatively large genome which were divergent to previously known viruses from the viral families Rhabdoviridae and Coronaviridae. Depending on specific characteristics of the target virus and the metagenomic community, different assembly and in silico gap closure strategies were successful in obtaining near complete viral genomes.

Viral metagenomic analysis of feces of wild small carnivores

Background: Recent studies have clearly demonstrated the enormous virus diversity that exists among wild animals. This exemplifies the required expansion of our knowledge of the virus diversity present in wildlife, as well as the potential transmission of these viruses to domestic animals or humans. Methods: In the present study we evaluated the viral diversity of fecal samples (n = 42) collected from 10 different species of wild small carnivores inhabiting the northern part of Spain using random PCR in combination with next-generation sequencing. Samples were collected from American mink (Neovison vison), European mink (Mustela lutreola), European polecat (Mustela putorius), European pine marten (Martes martes), stone marten (Martes foina), Eurasian otter (Lutra lutra) and Eurasian badger (Meles meles) of the family of Mustelidae; common genet (Genetta genetta) of the family of Viverridae; red fox (Vulpes vulpes) of the family of Canidae and European wild cat (Felis silvestris) of the family of Felidae. Results: A number of sequences of possible novel viruses or virus variants were detected, including a theilovirus, phleboviruses, an amdovirus, a kobuvirus and picobirnaviruses. Conclusions: Using random PCR in combination with next generation sequencing, sequences of various novel viruses or virus variants were detected in fecal samples collected from Spanish carnivores. Detected novel viruses highlight the viral diversity that is present in fecal material of wild carnivores.

Inferring Population Genetic Structure in Widely and Continuously Distributed Carnivores: The Stone Marten (Martes foina) as a Case Study

The stone marten is a widely distributed mustelid in the Palaearctic region that exhibits variable habitat preferences in different parts of its range. The species is a Holocene immigrant from southwest Asia which, according to fossil remains, followed the expansion of the Neolithic farming cultures into Europe and possibly colonized the Iberian Peninsula during the Early Neolithic (ca. 7,000 years BP). However, the population genetic structure and historical biogeography of this generalist carnivore remains essentially unknown. In this study we have combined mitochondrial DNA (mtDNA) sequencing (621 bp) and microsatellite genotyping (23 polymorphic markers) to infer the population genetic structure of the stone marten within the Iberian Peninsula. The mtDNA data revealed low haplotype and nucleotide diversities and a lack of phylogeographic structure, most likely due to a recent colonization of the Iberian Peninsula by a few mtDNA lineages during the Early Neolithic. The microsatellite data set was analysed with a) spatial and non-spatial Bayesian individual-based clustering (IBC) approaches (STRUCTURE, TESS, BAPS and GENELAND), and b) multivariate methods [discriminant analysis of principal components (DAPC) and spatial principal component analysis (sPCA)]. Additionally, because isolation by distance (IBD) is a common spatial genetic pattern in mobile and continuously distributed species and it may represent a challenge to the performance of the above methods, the microsatellite data set was tested for its presence. Overall, the genetic structure of the stone marten in the Iberian Peninsula was characterized by a NE-SW spatial pattern of IBD, and this may explain the observed disagreement between clustering solutions obtained by the different IBC methods. However, there was significant indication for contemporary genetic structuring, albeit weak, into at least three different subpopulations. The detected subdivision could be attributed to the influence of the rivers Ebro, Tagus and Guadiana, suggesting that main watercourses in the Iberian Peninsula may act as semi-permeable barriers to gene flow in stone martens. To our knowledge, this is the first phylogeographic and population genetic study of the species at a broad regional scale. We also wanted to make the case for the importance and benefits of using and comparing multiple different clustering and multivariate methods in spatial genetic analyses of mobile and continuously distributed species.

Amplicon metagenomics for elucidating microbial community changes associated to conventional and organic farming of Beta vulgaris and Solanum lycopersicum

Soil microbial community changes associated to conventional and organic farming of two relevant crops (Beta vulgaris and Solanum lycopersicum) were analysed through 16s rRNA amplicon sequencing. This study revealed microbial communities in the agricultural soils studied to be similar to other reported nutrient-rich microbiomes, and some significant differences between the microbial communities associated to the two farming practices were found. Some phyla (Chloroflexi and Thermi) were found to be present in different abundances according to soil treatment. As chloroplast interference can be a stumbling block in plant-associated 16s rRNA amplicon metagenomics analysis of aerial plant tissues, two protocols for bacterial cell detachment (orbital shaking and ultrasound treatment) and two protocols for microbial biomass recovery (centrifugation and filtration) were tested regarding their efficiency at excluding plant-DNA. An alternative method to the one proposed by Rastogi et al (2010) for evaluating the chloroplast-amplicon content in post-PCR samples was tested, and the method revealed that filtration was the most efficient protocol in minimising chloroplast interference.

Pseudohypoparathyroidism Type Ib Associated with Novel Duplications in the GNAS Locus

Context Pseudohypoparathyroidism type 1b (PHP-Ib) is characterized by renal resistance to PTH (and, sometimes, a mild resistance to TSH) and absence of any features of Albright's hereditary osteodystrophy. Patients with PHP-Ib suffer of defects in the methylation pattern of the complex GNAS locus. PHP-Ib can be either sporadic or inherited in an autosomal dominant pattern. Whereas familial PHP-Ib is well characterized at the molecular level, the genetic cause of sporadic PHP-Ib cases remains elusive, although some molecular mechanisms have been associated with this subtype. Objective The aim of the study was to investigate the molecular and imprinting defects in the GNAS locus in two unrelated patients with PHP-Ib. Design We have analyzed the GNAS locus by direct sequencing, Methylation-Specific Multiplex Ligation-dependent Probe Amplification, microsatellites, Quantitative Multiplex PCR of Short Fluorescent fragments and array-Comparative Genomic Hybridization studies in order to characterize two unrelated families with clinical features of PHP-Ib. Results We identified two duplications in the GNAS region in two patients with PHP-Ib: one of them, comprising similar to 320 kb, occurred 'de novo' in the patient, whereas the other one, of similar to 179 kb in length, was inherited from the maternal allele. In both cases, no other known genetic cause was observed. Conclusion In this article, we describe the to-our-knowledge biggest duplications reported so far in the GNAS region. Both are associated to PHP-Ib, one of them occurring 'de novo' and the other one being maternally inherited.

Identification of Sex and Female’s Reproductive Stage in Commercial Fish Species through the Quantification of Ribosomal Transcripts in Gonads

The estimation of maturity and sex of fish stocks in European waters is a requirement of the EU Data Collection Framework as part of the policy to improve fisheries management. On the other hand, research on fish biology is increasingly focused in molecular approaches, researchers needing correct identification of fish sex and reproductive stage without necessarily having in house the histological know-how necessary for the task. Taking advantage of the differential gene transcription occurring during fish sex differentiation and gametogenesis, the utility of 5S ribosomal RNA (5S rRNA) and General transcription factor IIIA (gtf3a) in the molecular identification of sex and gametogenic stage was tested in different economically-relevant fish species from the Bay of Biscay. Gonads of 9 fish species (, Atlantic, Atlantic-chub and horse mackerel, blue whiting, bogue, European anchovy, hake and pilchard and megrim), collected from local commercial fishing vessels were histologically sexed and 5S and 18S rRNA concentrations were quantified by capillary electrophoresis to calculate a 5S/18S rRNA index. Degenerate primers permitted cloning and sequencing of gtf3a fragments in 7 of the studied species. 5S rRNA and gtf3a transcript levels, together with 5S/18S rRNA index, distinguished clearly ovaries from testis in all of the studied species. The values were always higher in females than in males. 5S/18S rRNA index values in females were always highest when fish were captured in early phases of ovary development whilst, in later vitellogenic stages, the values decreased significantly. In megrim and European anchovy, where gonads in different oogenesis stages were obtained, the 5S/18S rRNA index identified clearly gametogenic stage. This approach, to the sexing and the quantitative non-subjective identification of the maturity stage of female fish, could have multiple applications in the study of fish stock dynamics, fish reproduction and fecundity and fish biology in general.